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The question of nuclear response functions in a homogeneous medium is examined. A general 
method for calculating response functions in the random phase approximation (RPA) with exchange 
is presented. The method is applicable for finite-range nuclear interactions. Examples are shown in 
the case of symmetric nuclear matter described by a Gogny interaction. It is found that the conver- 
. . . gence of the results with respect to the multipole truncation is quite fast. Various approximation 

' schemes such as the Landau approximation, or the Landau approximation for the exchange terms 

, only, are discussed in comparison with the exact results. 
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lO ' I. INTRODUCTION 

(N 

Infinite nuclear matter as a homogeneous medium made of interacting nucleons is not a system that can be ex- 
J> ' perimentally studied in the laboratory, but it is nevertheless a very useful and broadly used concept because of its 
\ relative simplicity and its connection with the inner part of atomic nuclei. In the remote environment of our planet 
' Earth this idealized system can modelize some parts of the compact stars. It is therefore important to have a micro- 
\ scopic description of nuclear matter based on nucleon-nucleon interactions. There are basically two main approaches, 
. either by starting from a bare two-body force and treating the many-body problem by Monte-Carlo methods or 
' Brueckner-Hartree-Fock method or using directly an effective nucleon-nucleon interaction adjusted to describe 

, the bulk properties of nuclear matter and finite nuclei in a mean field approximation. In the latter approach there are 
two types of interactions very widely used in a non-relativistic framework, namely the Skyrme-type forces 0, 01 and 
the Gogny-type forces . In this work we concentrate on the question of nuclear response functions using finite-range 
forces like the Gogny force. 

Q \ There are many physical issues that require the knowledge of the response function of the medium to an external 
^ . probe. Well-known examples are the electron scattering by nuclei or the propagation of neutrinos in nuclear matter. 
^ ■ In a mean field framework the response functions must take into account the effects of long-range correlations by the 
Random Phase Approximation (RPA) which is the small amplitude limit of a time-dependent mean field approach. 
For contact interactions of the Skyrme type the RPA response functions have been often studied (see, e.g., Ref. j3|)- 
On the other hand, RPA studies of nuclear matter with Gogny forces are more rare, and they usually involve some 
limiting assumption such as the Landau limit, or the small momentum transfer limit 

It is worthwhile at this point to clarify the terminology used in the literature. One often uses the short-hand name of 
RPA for the ring approximation of RPA which is obtained when the particle-hole (p-h) interaction is approximated by 
its Landau-Migdal form. Here, the main purpose is to treat exactly the exchange contributions of the p-h interaction 
and therefore, we keep the name RPA for this situation unless otherwise specified. This corresponds to what is called 
RPA with exchange (RPAE) in the electron gas physics. 

The method for solving the RPA equation with a finite range force is simple in principle. We show that, with 
a small number of terms in the multipole expansions the convergence is fast and the calculations are relatively 
easy. We also compare the exact RPA response functions with various approximations, namely the Landau-Migdal 
approximation made on the complete p-h interaction or on the exchange part of it. The latter approximation keeps 
the exact momentum transfer dependence of the direct interaction but it still has the simplicity of the Landau-Migdal 
treatment and therefore, it can be useful for very extensive studies of particle propagation inside matter. 

The outline of the paper is as follows. In Sec. II we present the general method for calculating RPA response 
functions with direct and exchange p-h interactions. In Sec. Ill we discuss the convergence of the multipole expansion 
using the Gogny force DIS. In Sec. IV we compare the exactly calculated RPA response functions with approximations 
of Landau-Migdal type. Concluding remarks are in Sec. V. 
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II. FORMALISM 

A general two-body interaction in momentum representation depends at most on 4 momenta. Because of momentum 
conservation there are actually 3 independent momenta. For the particle-hole (p-h) case we choose these independent 
variables to be the initial (final) momentum ki (k2) of the hole and the external momentum transfer q. This 
is illustrated by Fig. ^ We will denote by a = {S,T) the spin and isospin p-h channels with S=0 (1) for the 
non spin-flip (spin-flip) channel, and T=0 (1) the isoscalar (isovector) channel. The matrix element of the general 
antisymmetrized p-h interaction can be written as: 

(q + ki,kri|y|q + k2,k2-l) = ^l;(^")(q,kl,k2)p("^ (1) 

(a) 

where the projectors are p("'°' = l/g, P^^'°'> = (cti • crz)/^, P^"'^^ = (n • tz)/.? and P^^-^) (ai • a2){Ti ■ T2)/.g (the 
factor (7=4 is the spin-isospin degeneracy). 

q-l-K2 ^2 




q+ki kj 

FIG. 1: Direct and Exchange parts of the ph interaction. 

For a finite range interaction whose spin-isospin dependence is described by the usual Wigner, Bartlett, Heisenberg 
and Majorana terms, the components T^^^^ (q, ki, k2) are: 

V^'^\^MM) = i?("'/(9)+i?^"^/»(ki-k2) + i?(") . (2) 

The direct terms D^"^f(q) depend only on the modulus q while the exchange terms p("^/i(ki — k2) depend on ki — k2. 
The last term R^"^ accounts for rearrangement contributions and it must be included if the starting nucleon-nucleon 
effective interaction has a density dependence (Tll |. 

As an example, let us consider as the starting interaction V the effective Gogny force @ which is often used for 
nuclear matter and nuclear structure studies. It consists of a sum of two Gaussians having different ranges and spin- 
isospin dependences supplemented by a contact term depending on the local density. For such a force, the functions 
f{q) and h(ki — k2) are: 

fiq) = e-y^^' 

/i(ki-k2) = e-^f^^^-"-)''^' , (3) 

where ^ is the range parameter of the Gaussian form factor. For a Gogny-type force, the direct and exchange 
contributions to Eq. (PJ are obtained by summing over the two Gaussians. The expressions of and P'") in terms 
of the spin-isospin coefficients and the range of the Gaussian, as well as the rearrangement term P(") are shown in 
Table 1. 





(0,0) 


(0,1) 


(1,0) 


(1,1) 


^(S,T) 


7rV^^^(4W + 2B-2H - M) 
■Ky/^^?{-W -2B + 2H + 4M) 
i(7 + l)(7 + 2)t3P^ 


TT^fi'\-2H - M) 

-Ky^fl\-W -2B) 
^{l + 2X3)t3p'' 


TT^fi^{2B - M) 
-K^^i^-W + 2H) 

- (1 - 2X3)t3p'' 


■K^/^^JL^{-M) 
TV^lJ^i-W) 

-tap'' 



TABLE I: Direct, exchange and density-dependent components of the p-h interaction corresponding to a Gogny-type interaction. 
W,B,H,M,p,t3,j are interaction parameters, p is the local density. 

Let us consider for simplicity an infinite nuclear medium at zero temperature and unpolarized both in spin and 
isospin spaces. At mean field level this system is described as an ensemble of independent nucleons moving in a self- 
consistent mean field generated by the starting effective interaction treated in the Hartree-Fock (HF) approximation. 
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The momentum dependent HF mean field, or self-energy determines the single-particle spectrum e(fc) and the Fermi 
level e{kp). 

To calculate the response of the medium to an external field it is convenient to introduce the p-h Green's function, 
or retarded p-h propagator G(q, a;, ki). From now on we choose the z axis along the direction of q. In the HF 
approximation, the p-h Green's function is the free retarded p-h propagator [l^ : 

w + e(fci) - e(|q + ki)| + t?/ 

It is customary to go beyond the HF mean field approximation and to take into account the long-range type of 
correlations by resumming a class of p-h diagrams. One thus obtains the well-known Random Phase Approximation 
(RPA) whose correlated Green's function G((3', w, ki) satisfies the Bethe-Salpeter equation: 

■^Vph{qMM)G{q,u;M) ■ (5) 
Finally, the response function xil,^) in the infinite medium is related to the p-h Green's function by: 

j^^G{q,uj,'ki) , (6) 

Equations should be understood for each {S, T) channel. The Lindhard function x^"^ is obtained when the free 

p-h propagator G^"^ is used in Eq. J^J. 

If the p-h interaction Eq. (j5Jl is treated in some approximation so as to simplify its ki — k2 dependence, then 
solving Eq. lO can be made easier. For instance, neglecting the exchange terms or treating them in Landau 

approximation and keeping only I = terms lead to the familiar ring approximation of RPA, and the response 
function © can be expressed in terms of the Lindhard function x'°' . In the past there have been studies of nuclear 
matter response functions with Gogny-type interactions under simplifying assumptions. For instance, RPA response 
functions in the long-wavelength limit, i.e., for vanishing momentum transfer q have been calculated in Ref. 
Actually, that study was done in the Landau approximation whereas we are looking for a complete RPA calculation 
with finite range interactions. 

The method presented in this paper allows one to obtain RPA Green's functions for all values of q and lo without 
further approximation. This will be useful for studying processes such as the propagation of particles inside nuclear 
matter. It is a cumbersome task to solve directly in the 3-dimensional momentum space the Bethe-Salpeter equation 
with the full p-h residual interaction. The general method of solution proposed here is to expand the Green's functions 
and the p-h interaction on a complete basis of spherical harmonics and to transform Eq. © into a set of coupled 
integral equations on the radial momentum (i.e., the momentum modulus) variable. The expansion of the unperturbed 
p-h Green's function is simple because it has no dependence on the angle (j) of the vector ki: 



G("'('7,'^,ki)=^Gf (<?,^,fci)rLo(fci) . (7) 



^^(0) I 
L 

As for the residual interaction Eq. (|2j it depends on the modulus |ki — k2| and therefore, its expansion is: 

Vj,h{qMM) ^Y.^L{q,h,k2)Yl,,{h)YLM{k2) ■ (8) 
LM 

Because of the structure of the multipole expansion of Vph, the expansion of the RPA Green's function is similar to 
that of G^o): 

G(g,a;,ki) = ^GL((7,c^,fei)rLo(fci) • (9) 

L 

Indeed, by inserting back this expansion of the propagator into Eq. lO, a consistent result is obtained. Making use of 
the multipole expansions II7I9II one can transform the Bethe-Salpeter equation into a set of coupled integral equations 
for the multipole components of the RPA Green's function: 

GL{q,u;,k,) = G'-°\q,LJ,ki)+J2 J '^^FLL'{q,u;;kuk2)GL'{q,0J,k2) , (10) 
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where we have defined 



where (LiMiL2-^2|-£'3-^3) is a Clebsh-Gordan coefficient. The angular momenta L and L' entering the above expres- 
sions are unlimited in principle. We will see in the next section that in practice the convergence is very fast and very 
few terms are necessary to obtain a good accuracy. Finally, the response function Eq. J^J can be expressed as: 



\/4tt f 

Xiq,^) = J^^S J GL=n{q,w,k)k'^dk , (12) 

where only the L ~ multipole of the RPA Green's function is required. However, one has to solve the full system of 
coupled equations Ijl0|l . since the interaction couples diflerent multipoles. In practice, a complete calculation implies 
the choice of a cut-off value Lmax for the summations on angular momenta, and a grid of points in momentum space 
in order to transform the integrals into discrete sums. Then, Eq. (fTTHl is solved by a matrix inversion. For the results 
presented in the next sections we have chosen a grid with a constant number of points. With 100 points, we have 
obtained a good comparison of the free response function with its analytic form, the Lindhard function. The upper 
limit is kmax — kp + q. Hence, for q = fcF/10, Afc=0.017 fm"^ and kmax='^-T^ fm^^ and for q = kp, Afc=0.031 fm~^ 
and fcma2;=3.14 fm"^ 



III. SYMMETRIC NUCLEAR MATTER RESULTS 



The Gogny force DIS has been chosen to discuss some results calculated in symmetric nuclear matter. First, the 
unperturbed p-h Green's function G^^^ (g, k) is calculated using the Hartree-Fock solution corresponding to DIS. 

It is just the familiar Lindhard function and the multipoles Gf\q,u;,k) can be easily calculated numerically. Then, 
an important practical issue is the sensitivity of the solution of Eq. Ill()|l to the number of multipoles included in the 
calculation. To examine this point we have performed calculations with different values of Lmax- In Figs. I2I3I the real 
and imaginary parts of the RPA response function (q, oj) are displayed for several values of L^ax and two values 
of q, namely 27 MeV (~ kp/lQ) and 270 MeV (~ kp). 




20 -20 

0) [McV] CO [MeV] 



FIG. 2: Real and imaginary parts of the response function x for q=27 MeV. Dotted: Lmax = 0. Dashed: Lmax ~ 1. Solid: 
Lm.ax = 2 (converged). 



It can be seen that for small q the convergence is reached for Lmax='^ (the solid and dashed curves are indistin- 
guishable in the figure). Increasing the value of q, one sees that there is still a small difference between Lmax = i- and 
2. The convergence is reached in fact for L,nax=3. 

A quantitative measure of the degree of convergence can be provided by the expected symmetries of the real and 
imaginary parts of the response functions. Indeed, the former should be symmetric and the latter anti-symmetric 
with respect to ll)=0. Let us introduce for this purpose the following symmetry parameters: 



/ du;\Rcx{uj)~Rcx{~u;)\/ f , 
Juj>Q I JLu>a 



du} |Rex('^) + Rex(— w) 



/ dbj \lmx{uj) + lmxi^'^)\ / / dw |Inix(w) — Inix(— Ct;)| 

Juj>0 I Juj>0 



(13) 
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FIG. 3: Same as Fig.l for g=270 MeV. 

The results corresponding to the cases shown in Figs. H2I8II are indicated in Table lUl Although the mentioned 
symmetric behavior is seen in the figures, the results of TablelTTI allow one to conclude more quantitatively that, in the 
range of momentum transfer up to kp, the symmetry criteria are satisfied within 1% using Lmax=S. In this case, the 
matrices involved in the solution of Eq. (jlOII have relatively moderate sizes (less than 500x500) and the calculations 
are rather fast. A similar numerical study for the response functions in the channels other than {S = 0,T ~ 0) has 
been performed, leading to the same type of convergence as a function of Lmax- 







Lmax 




q (MeV) 





1 


2 


3 


4 


5 




27 


2.10 


2.27 


0.42 


0.06 








270 


14.7 


14.0 


3.75 


0.94 


0.56 


0.52 




27 


3.60 


2.32 


0.60 


0.12 








270 


14.0 


13.0 


3.38 


0.85 


0.36 


0.35 



TABLE II: Symmetry parameters corresponding to the cases shown in Figs. I2I3I 

To understand this rapid convergence, it is necessary to analyze in some detail the multipoles of the p-h interaction. 
To this end, we plot the dimensionless multipoles 

or _L 1 

VL{q,h,k2) = — NoVL{q,kuk2) (14) 

for a fixed value of q in the plane {ki/kF,k2/kF) of hole momenta. In this expression, Nq = gmfkr / ['^t^'^) is the 
density of states, and the factor multiplying Vl fixes the scale: the value at (7=0, fci = fc2 = kp gives the familiar 
Landau dimensionless parameter. The factor 47r comes from the use of spherical harmonics Ylo instead of Legendre 
polynomials Pl in the multipole expansions. 




kj/kp kj/kp 



FIG. 4: The dimensionless monopole VQ^''^\q, ki, as a function of the momenta k-i/kp and k'z/kp at p = po. The left panel 
represents the case g=0 and the right one g=270 MeV. 
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The monopole L=0 case is plotted in Fig.^lfor q=0 (left panel) and (7=270 MeV (right panel). The effect of a finite 
transferred momentum q only affects the monopole component of the p-h interaction (see Eq. [JJ , and it produces an 
overall translation of the interactions drawn in the figure. For this specific channel, it induces a repulsion. 




kj/kp kj/kp 



FIG. 5: The dimensionless multipoles Vf^ ' (q,fci,fc2) as a function of the momenta ki/kp and k2/kF at p — po- The left panel 
represents the case L=l and the right one L=2. 

The multipoles L=l and L=2 are plotted in Fig. [3 in the left and right panels respectively. This series of figures 
for L=0,l,2 shows explicitly that the multipole expansion of the p-h interaction is rapidly converging after the first 
few angular momenta. 

Inside the {ki,k2) domain considered, the variations of the multipoles Vl span the range of values (-2:1) for L=0, 
(-1:0) for i=l, (-0.5:0) for L=2, (-0.3:0) for L=3 while beyond L=4, the multipoles are practically negligible. In 
conclusion, the full convergence is achieved using Lmax=S. 

IV. LANDAU APPROXIMATIONS 

Once the convergence of the method has been proved, it is useful to analyze some approximations employed to 
obtain the response function. The Landau-Migdal approximation is often used in the literature because it 
simplifies greatly the calculation of RPA response functions. This approximation was used in Ref. for the Gogny 
interaction Dl, based on the solution of the kinetic equations. We should mention that no analysis of convergence 
was made. The approximation consists in assuming that the interacting particle and hole are on the Fermi surface 
and that the interaction takes place only in the limit q = 0. That is to say that each multipole of the p-h interaction 
is replaced by the constant Vl{0, kp, kp). The validity of this assumption can be checked by inspecting Figs. 14151 For 
the considered interaction and (5, T) channel, the first Landau dimensionless parameters are Fo=-0.38, Fi=-0.91 and 
F2=-0.33. It can be seen in these figures that the corresponding multipoles of the p-h interaction are far from taking 
a constant value. It is well known that the validity of the approximation is limited to very small values of g, because 
in this case the physically relevant values of fci and k2 remain around the Fermi momentum. However, it is worth 
keeping in mind that in what concerns the response function the differences between the true and the approximated 
p-h interactions are smeared out by the kinetic constraints of the phase space. 

The g=0 hypothesis can be easily relaxed. Indeed, Eq. ^ shows that only the exchange of the p-h interaction 
depends on ki and k2- Thus, we could keep unchanged the direct term D'^'^^f^q) of the interaction Eq. ® with its 
full g-dependence and make the Landau approximation on the exchange term i?("'ft,(ki — k2) only. We denote this 
choice as the Landau Approximation For Exchange Term (LAFET). One may hope that this procedure will improve 
the usual Landau approximation by treating approximately the exchange term only. The simplicity of the Landau 
approximation is preserved, the only change being that the monopole contribution acquires a g-dependence coming 
from the direct term. 

The general method presented in Sec. ^ can be easily applied to both approximations. In the appendix El it 
is shown that if the p-h interaction is independent of the hole momenta fci and k^, the system of coupled integral 
equations Eq. ifTTijl for the multipoles Giiq, i^, k) can be transformed into a set of algebraic equations for their integrals 
Jdkk^GLiq,uj,k). 

We represent on Figs. 16171 a comparison of the structure function 

^(^.i^)(g,a.) = --W^^^)(g,t^) (15) 
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FIG. 6: Structure function for g=27 MeV and L^ax = 3. Solid thin curve: HF. Dotted: Landau. Dashed: LAFET. 

Solid thicli line: converged. 




GO [MeV] GO [MeV] 

FIG. 7: Structure function S'^'^' (g, a;) for g=270 MeV and Lmax = 3. Solid thin curve: HF. Dotted: Landau. Dashed: 
LAFET. Solid thicli line: converged. 

extracted from a converged solution (solid thick curve) with Landau (dotted) and LAFET (dashed) approximations 
for q = kp /lO and q = kp. The HF solution (solid thin curve) is also displayed as a reference. All the curves shown 
correspond to L,nax=3, to guarantee the convergence of the solution as shown in Sec. IIIII As expected, the Landau 
approximation is a good one for small transferred momentum q and all {S, T) channels. However, for the highest 
transferred momentum considered in this paper, the validity of the Landau approximation and LAFET response 
functions is very {S,T) dependent. As a rule, LAFET induces a response closer to the exact one. Hence, LAFET 
could be considered as a very simple extension of the Landau approximation which allows one to evaluate processes 
which involve finite transferred momenta. 

There are situations where it is necessary to calculate accurately response functions of a nuclear medium to an 
external probe, as we already mentioned in the introduction. When the transferred momenta are not small compared 
to the Fermi momentum, one must perform the full calculations with the method presented in Sec. II. Alternatively, 
it is possible to use the LAFET method, which is very simple and efficient if extensive calculations are needed, and 
improves the Landau approximation. 

V. CONCLUSION 

The main purpose of this work is to present a general method for obtaining nuclear response functions in an infinite 
medium within a Hartree-Fock-RPA framework. Starting with finite range effective interactions like the Gogny 
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interaction, approximate methods using the Landau approximation are available in the literature, but surprisingly 
not much beyond this approximation can be found. The method proposed here simply consists in expanding the Bethe- 
Salpeter equation onto a spherical harmonics basis and therefore, the calculations can be carried out in principle up 
to any degree of accuracy if one includes a sufHcient number of partial waves. In practice, the case study that we 
have discussed in this work shows that the convergence is very fast and that the number of multipoles needed is very 
small. This result holds not only for small values of momentum transfer but even at values in the range of the Fermi 
momentum. The fast convergence is related to the properties of the effective p-h interaction. 

This general method of solving the Bethe-Salpeter equation suggests also an approximation scheme beyond the 
standard Landau approximation, the LAFET scheme where the full g-dependence is kept in the direct p-h interaction 
and the Landau approximation is done only on the exchange p-h interaction. The standard Landau approximation and 
LAFET are compared with the exact response functions and it is shown that the LAFET results show an improved 
agreement with the exact results. This approximation can be useful for extensive calculations when the numerical 
effort required by exact calculations becomes heavy. 
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APPENDIX A: LANDAU APPROXIMATIONS 



Let us show that if the p-h interaction is independent of the hole momenta ki and k2 , the system of coupled integral 
equations Eq. for the multipoles Gl (g, w, fc) can be transformed into a set of algebraic equations for the following 
quantities 



RLiq,uj)=g^ J dkeGL{q,LO,k) . (Al) 

The factors before the integral are chosen such that the response function is given by Rl=o- 

Let us consider the multipoles of the p-h interaction, Eq. in the particular case where the hole momenta ki 
and k2 lie on the Fermi surface. For this specific interaction, the functions Fll' entering Eq. no longer depend 
on k2- Since all the momentum dependence is now contained in the p-h propagators, Gl and G^^\ integrating over 
ki the system Eq. l(T(1| can be transformed into the set of equations: 

Rl [q, to) = i?f iq.io)+J2 All' (<Z, uj)Rl' (g, u;) , (A2) 

LU 

where Rf/' is defined as in Eq. (IA1|I for the free propagator and 

ALL'iqM = ^- L[ (2L^ + i)(2A+l) j (^0^01^0) R, . (A3) 

It is worth noting that the Landau parameters are given by Jl = VL{q = 0, kp, kp)/No. 

As an example, we give here the response function in the Landau approximation with Lrnax='^'- 

AO) 

where 



,2 >- ' y s '10^ 

iFi)(l + i[-9i. , 5 

plays the role of the induced p-h interaction. In this expression, we have defined ly ~ ujm* / (qkp), and Fi = fiNo are 



ti inter; 

the dimensionless Landau parameters. 
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